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SUMMARY 


The paper treats the general problem of solving for the steady state (time 
domain) dynamic response (i.e., NASTRAN rigid format-8) of a general elastic 
periodic structure subject to a phase difference loading of the type encountered 
in traveling wave propagation problems. Two types of structural configurations 
are considered; in the first type, the structure has a repeating pattern over a 
span that is long enough to be considered, for all practical purposes, as 
infinite; in the second type, the structure has structural rotational symmetry 
in the circumferential direction. Due to the periodic nature of the structure 
and the traveling wave characteristics of the loading, one need only "cut out" 
and subsequently model a typical periodic region of the total structure, 
wherein appropriate periodic boundary conditions (i.e., unknown forces and dis- 
placements are forced equal, except for unknown phase angle, for corresponding 
points on both cuts) are used along the cuts. The paper presents both the 
theory and a corresponding set of DMAP instructions which permits the NASTRAN 
user to automatically alter the rigid format-8 sequence to solve the intended 
class of problems. The new input to a standard version NASTRAN run is a set 
of alter cards, PARAM cards, and direct input matrix (DMI) partitioning arrays 
which are needed for the purpose of partitioning and correspondingly restruc- 
turing the internal NASTRAN mass, damping and stiffness matrices. Final 
results are recovered as with any ordinary rigid format-8 solution, except that 
the results are only printed for the typical periodic segment of the structure. 
A simple demonstration problem having a known exact solution is used to illus- 
trate the implementation of the procedure. 



SYMBOLS 


[B] 

Damping matrix of n in periodic substructure 

{F} 

Total applied force vector 


i 

V ' 1 


[I U ] 

Diagonal unit identity matrix 


I 

Number of degrees-of-freedom for 

interior nodes 

[K] 

j. u 

Stiffness matrix of n tn periodic 

substructure 
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L Spatial period of substructure 

r 

L Number of degrees-of- freedom for left cut nodes 

1. L 

[M] Mass matrix of n tr periodic substructure 

R Number of degrees-of-freedom for right cut 

t Time 

{ U> Displacement vector 

x Spatial coordinate 

co Angular driving frequency 

0, ip Angle of incident wave 

y Phase constant (complex form) 

y* Phase constant ( real part) 

INTRODUCTION 


A periodic structure consists of a number of identical substructures, 
coupled together in identical manners to form the whole system, see for example 
figure la. For such systems, under certain loading conditions, it is often 
possible to treat only one representative substructure in order to obtain the 
general response for the whole system. For example, if the loading is exactly 
the same for all substructures, the latest (and even some earlier) versions of 
NASTRAN can directly solve this class of problem for both static and steadystate 
cases ( i . e . , rigid formats 1 and 8). In the case of steady state dynamics 
problems (rigid format-8) involving traveling wave propagation type inputs, 
there is a slightly more general loading condition on each periodic substruc- 
ture, namely that the loading on each substructure is identical except for a 
known phase constant y . More specifically, the relation between the applied 
force vector {7} n in the n th substructure and the one, {7} n+ j in the n™ sub- 
structure is given by 


{T} 


n+1 




(1) 


where y is a known phase constant. For the class of problems addressed in this 
paper, the phase constant is a purely imaginary constant, i.e». 


y - 0.0 + iy* 


( 2 ) 



and physically refers to the fact that there is no difference in energy loss in 
processing results from one substructure to the next. 

Brillouin (ref. 1) points out that wave motion in periodic systems have 
been studied for nearly 300 years, wherein physicists and electrical engineers 
have worked in this field in problem areas relating to optics, crystals, elec- 
trical transmission liner, etc. (Elachi, ref. 2, provides a comprehensive list 
of 287 references in this field). Applications of this theory to engineering 
structural analysis and solid mechanics type problems is only recent. Refer- 
ences 3, 4, and 5 are typical of analytical solutions to this type problem for 
simple configurations consisting of beams, grillages and plate structures. 
References 6 and 7 represent a significantly more general approach to the prob- 
lem wherein their application of the theory of finite elements enables one to 
solve a much larger class of problems involving rather arbitrary structures 
than one could treat by purely analytical techniques. References 6 and 7 
appear to restrict themselves to the problem of determining the conditions 
(i.e., values of the steady state response freqency w) under which propagating 
or non-propagating free wave motion will occur in the absence of explicit 
external driving forces. 

In the work presented here, periodic structures with explicit external 
driving forces satisfying Equations (1) and (2) are applied to each substructure 
as illustrated in Figure la. 

If the loading and spatial boundary conditions on each substructure are 
the same, except for the phase difference, y*, (i.e., the loading for two 
typical points in two neighboring substructures separated by the period length 
Lp,are the same except for a multiplying factor of e 11J *) it follows that the 
response in each substructure is also the same except for the phase difference 
y*. A simple example of such a case is a propagating pressure wave passing 
across an infinitely long ribbed plate as illustrated in figure 2b (the plate is 
in air and no air-structure interaction effects are considered). The propagat- 
ing surface loading wave is given by the formula 


p » p^'fkx + ut) 


(3) 


thus the phase difference between any two neighboring substructures is y* = kLn, 
where Lpis the spatial period of the periodic system, p Q is the input loading 

pressure amplitude, to is the steady state driving frequency, x is the horizontal 
spatial coordinate, and k is the wave number of the loading wave. 

Other examples of the phase constant relative to a particular example are 
shown in figure 2. In figure 2a, we have a known pressure wave loading propa- 
gating parallel with the axis of the ribbed cylindrical shell; here, the phase 
constant y* is analogous to the figure 2b example and needs no further explana- 
tion. In figure 2c, the incident wave is incident at an oblique angle 0 and 
is incorporated into the formula for the phase constant given in the figure lc. 


133 



There is a special case where the ends of the structure do not extend to 
infinity (i.e., the ends never meet) but instead are connected cyclically, as 
in figure 3 for example. For such cyclical cases, y* must satisfy an addi- 
tional constraint, namely y* = 2-fr/n where n = 1, 2, 3,* * * . For example, in 
the figure 3 case, n = 8, thus y* = tt/4. 

We limit ourselves to problems having "one-dimensional periodicity", 
whereby this term we imply that only two cuts are needed (we shall refer to 
these as the left and right cuts, see figure lb) to separate the typical sub- 
structure from the system. The response within such a substructure can be, 
however, mul ti -dimensional . The remainder of the paper focuses on the proce- 
dure for obtaining the displacement and stress response within one typical 
block of the periodic system. The typical substructure block can be made up of 
various types of structural elements (including both elements with structural 
damping and nodes with scalar dampers attached) contained within the NASTRAN 
library of elements (e.g., CQDMEM, CQDMEM1, CBAR, CONROD ... etc.). 


SOLUTION FORMULATION 


The solution procedure presented here is very similar to the one in refer- 
ence (6), except for the fact that here we are considering problems with expli- 
cit forcing functions. The first step in the solution procedure is to "cut 
out" the typical substructure from the overall periodic structure as illus- 
trated in figure lb and to subsequently replace the cut nodes with the internal 

forces ( {"F? } , for the left cut and {F^} for the right cut) that existed at 

X/ l 1 ill 

those nodes before cutting. The displacements at the cut nodes are similarly 
denoted by {U 0 } and {U } where subscripts a and r denote left and right and 

the subscript n denotes the n tn substructure. Since we are only focusing on 

the results for the n tn substructure, it is convenient to drop the subscript 
n from here on for notational convenience. 

The governing equations of motion for the substructure are first expressed 
in the familiar finite element form 


[M]{U> + (BHU) + (KHU) + {F} (4) 

i. i- 

where [M], [B], [K] are the mass, damping and stiffness matrices of the n tr 
substructure respectively; {U} is the displacement vector of all nodes of the 

j_ i- 

n u substructure; {F} is the generalized force fector; (*) = d( )/dt, and bars 
above the variable denote the fact that the variables are complex and that the 

harmonic time response e lwt has not yet been suppressed thus 

{IT} = {U}e lwt {F} = (F}e iwt (5) 


134 



The next major step is to partition the matrices and vectors of equation 
(4) into left cut unknowns, right cut unknowns and interior unknowns (subscripts 
&, r, i refer to left, right and interior respectively and L, R, I refer to the 
total number of displacement component unknowns for the left, right and interior 
domain respectively; note due to periodicity, L = R. Thus it follows that after 
partitioning we have 


[M] = 


M 


in 




W 


( M £i I *tr 
M ii I M ir 

■ - H - - 

M . I M 
n 1 rr 


(L+R+I) x (L+R+I ) 




[B] = 


B u 


I B£i _[ B ^ r 

B ii I B ir 


rl 


B . I B 
r i 1 rr 



!< U 


K 0 “ 

Jlr 

[K] = 

K i* 

K 1i 

| M — M 

K ir 


Kr* 

l K H 

1 

K rr 




_ 



(1) x (L+R+I) 


( 6 ) 



Note the generalized force vectors {F} has been further decomposed as the sum of 

an unknown force vector, {T 0 }, (which denotes the yet unknown internal forces 

existing at the cuts in the structure) and a known applied force vector, {T 9 }, 
(which denotes all known forces existing within and at the cuts of the periodic 
substructure). 
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The full periodic structure is cut (left and right cuts), therefore it 
follows that the internal nodal forces normally existing at the cuts now play 
the role of external (as yet unknown) applied forces. 

The special case of an externally applied force appearing at a left or 
right cut requires special attention in that one component of the total force 
vector is due to the externally applied force and the other component is due to 
the internal force at the cut. The external force value on a cut must be 
shared between the generic substructure block being analyzed and its immediate 
neighbor; consequently these end type external force values are divided in half 
(see, for example, the situation in figure lc where node 1 lies on the left 
cut). 


A further relation that is needed in the formulation relates to the fact 

J. L 

that the right end of the n tn generic substructure is the beginning (left end) 
of the n+1 generic substructure, thus from Equation (1) it follows that 


{1 Y } = " eU{ ^ } 

<V - eU <V 


( 7 ) 


where the minus sign in the first of Equation (7) accounts for the fact that 
internal nodal forces acting as external forces on the right cut of generic 
substructure n are opposite in sign to the internal nodal forces acting on the 
left cut neighboring substructure n+1. 

The next step in the development is to substitute Equations (5) and (6) 

into Equation (4); the subsequent cancellation of e 1a)t permits us to drop the 
bar superscript notation thus arriving at a "reduced form of Equation (4)". 

At this point there are five groups of unknowns, namely {U„> , {ll 0 }, {U.}, {F^>, 

{F^}. The three row partitions of the reduced Equation (4) in conjunction with 
the two Equations (7), provide 3+2=5 corresponding groups of equations to 
balance the five groups of unknowns. Next, we substitute Equations (7) into 
the reduced form of Equation (4), and subsequently employ the third row parti- 
tion of reduced Equation (4) to eliminate the {F^} unknown. Doing these opera- 
tions result in the following set simultaneous equations for the displacement 
unknowns {U^} and { } : 
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-0) 2 

[MLL] 

+ 1 

-W 2 

[ML I] 


ico 

[BLL] 

+ 1 

iw 

[BLI] 

+ 

i— — 

[KLL] 

1 

- 1 

wm — 

[KLI] 

.. 

-U) 2 

[MIL] 

+ | 

-w 2 

[Mil] 

+ 

iw 

[BIL] 

+ 1 

iw 

[B 1 1 ] 

+ 


[KIL] 

1 


[KII] 

- 


(L+I) x (L+I) matrix 


tuy 


{fJ} + e" y * i {F^} 

< 


. ------- ► 

<V 


{F?} 



- 


(1) x (L+I) vector 


( 8 ) 


where 


[MLL] = [M^J + Cosu* [M £r J + Cosy* [M^] + [M rf ] 

[BLL] = [B u ] + Cosy* [B £f ] - a>Siny* [M^] + (Si ny*/w) [K^] 

+ Cosy* [B rA ] + [B rr ] + wSiny* [M^] - (Siny*/w) [K r£ ] 

[KLL] = [K u ] + Cosy* [K^] - wSiny* [B^] 

+ Cos w* tK r y + [K rr ] + wSiny* [B^] 

[MU] = [M Ai ] + Cosy* [M H ] 

[BLI] = [B^] + Cosy* [B ri ] + wSiny* [M ri ] - (Siny*/w) [K H ] 

[KLI] = [K £ .] + Cosy* [K f1 ] + wSiny* [B^] (9) 

[MIL] = [M u ] + Cosy* [M- r ] 

[BIL] = [B u ] + Cosy* [B- r ] - uSiny* [M. r ] + (Siny*/ W ) [K^] 

[KIL] = [K^] + Cosy* [K ir ] - wSiny* [B ir ] 

[Mil] = [M..] 

[B 1 1 ] = [B..] 

[KII] = [K^.] 

At this point, the linear set of complex algebraic Equations (8) can be solved 
for the unknown displacements {U^}, (U.. }. The unknown displacement at the 

right cut, {U r > can be easily computed with the second of equation (7). The 
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size of the algebraic system is governed by the (L+I) x (L+I) coefficient 
matrix (i.e., matrix multiplying the unknown displacement vector) where L+I 
equals the number of left cut, L, plus interior, I, unknown displacement com- 
ponents. Typically I»L, therefore it is not a very big additional burden on 
the equation solver to include the second of Equation (7) as part of the overall 
system (actually, we add R extra unknowns, {U r > , and R extra equations (where 

R = number of right cut unknowns). Thus in place of Equation (8), we consider 
the slightly larger, but equivalent system of 


-w 2 [MLL] + 
iw [BLL] + 
[KLL] 


_ _ _ j _ _ 


-w 2 [MIL] + 
iw [BIL] + 
[KIL] 


_ j 


iw [BRL] + 
[KRL], 


-w 2 [ML I] + 
iw [BLI] + 
[KLI ] 


-w 2 [Mil] + 
iw [BII] + 
[KII] 

[ 0 ] 


[ 0 ] 


[ 0 ] 


j _ _ 


[KRR] 


(L+I+R) x (L+I+R) 
coefficient matrix 


<v 


<F*> ♦ 


e-^F*} 



r no) 


where [BRL] = (siny*/w) [I u ] 
[KRL] = Cosy* [I u ] 
[KRR] = - [I U ] 


( 11 ) 


and 


[ I u ] = diagonal unit identity matrix 


The RI block of the displacement coefficient matrix in Equation (10) above 
is identically zero, thus the bottom R rows of the system of simultaneous equa- 
tions are totally independent of the solution to the top L+I rows. The length 
of the solution vector L+I+R is of exactly the same length of the original sub- 
structure matrix (Equations (4) and (6)), consequently the modification of the 
DMAP instructions becomes simplier because of the fact that one need only 
intercept the logic of the equation solver and replace the existing mass stiff- 
ness and damping matrices with the modified matrices defined by the new coef- 
fient matrix of Equation (10) (and associated new entry definitions from Equa- 
tions (9) and (11)). Since the length of the solution vector is still the same 
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as the original problem before modification, the post processing DMAP opera- 
tions for displacement printout, stress recovery, etc. need not be modified. 

An alternate scheme (although not yet implemented) would be to modify the 
input to the complex equation solver to accept the smaller (L+I) x (L+I) coef- 
fient matrix used in Equation (8) directly. After solving the smaller L+I 
length displacement vector, the full vector (i 0 e«, attaching the missing {U^} 

portion) can be formed by expanding it to length L+I+R via the second of Equa- 
tion (7). Finally, stress and displacement results can be processed in the 
usual way with existing DMAP operations. 


RIGID FORMAT -8 DMAP MODIFICATION FOR NASTRAN 


The periodic structure capability described in the previous section can be 
implemented in a standard version of NASTRAN. In particular, the DMAP sequence 
required to perform the necessary operations are listed in Appendix A. This 
DMAP sequence was checked out on an 1108 computer, standard version of level 
15.5 NASTRAN and is introduced in the EXECUTIVE CONTROL deck with the following 
instructions: 

ALTER 138 

(see Appendix A for specific instructions) 

ALTER 139,139 

(replace KDD, BDD, MDD with KDDX, BDDX, MDDX 
within call arguments of FFRD module level 15.5 

see Appendix A for detailed instruction) implementation 

ALTER 140 

(Conditional print statement, see Appendix- 
A for detailed DMAP instructions) 

These same level 15.5 DMAP instructions can also be applied to level 17.0 
NASTRAN, the only difference being that 

replace ALTER 138 with ALTER 158 

replace ALTER 139,139 with ALTER 159,159 

replace ALTER 140 with ALTER 160 

It is pointed out, however, that the level 17.0 modifications described 
above have not actually been tried although due to the similarity of the change, 
the DMAP sequence is expected to work. 
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It is important to note that we are modifying the standard NASTRAN unknown 
displacement vector coefficient matrix just prior to the entry into the FFRD 
module used for the solution to the simultaneous complex algebraic equations. 
The implication of this statement is that the row numbering scheme for the dis- 
placement vector has already accounted for the fact that single point con- 
straints, multipoint constraints and omitted coordinates have already been 
accounted for. Thus, for example, the length of the {U^} vector, L, is not 

simply the number of nodes on the left cut times the degrees-of-freedom oer 
node, but rather is less by the amount corresponding to the number of SPC's, 
MPC‘s and OMIT's relating to the nodes along the left cut. Similar comments 
apply to the length of the {U^} and { U r } vectors. The understanding of the 

above displacement vector length comments must be clearly understood by the 
user before attempting to fill out the input data matrix partitioning vectors 
CV100, CV010, CV001 defined later in this paper. 


INPUT DATA FOR NASTRAN RUN 


The BULK DATA input to a typical periodic structure run consist of two 
basic parts. The first part corresponds to the usual bulk data input cards 
normally required to make a NASTRAN run (e.g., GRID CARDS, ELEMENT CARDS, DAREA 
CARDS, FREQ CARD, DLOAD CARD, etc.); the second part consists of special input 
cards that are explained in the following text. 


PARAM Cards 


These cards are used to enter various matrix coefficients appearing in 
Equations (9) and (11); especially the constants 


Text Variable 

Computer Variable 

Cosy* 

s CMST 

0.0 + i co 

5 FIj*foEG 

- 03 2 

= Nj2fMEG2 

-1.0 

= Nj^NE 

-siny*/w 

= NSMSBpf 

-(siny*)to 

= NSMST^ 

+1.0 

5 ^ note 

s i ny*/w 

= SMSBJ^ 

(siny*)w 

= SMSTtf 


0 = zero 
0 = letter 


(12) 
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are read in on standard NASTRAN PARAM cards where y* is the phase angle defined 
in Equation (2) and u is the angular driving frequency in radians per second 
(00 = f • 2 tt where f = driving frequency in H z specified on the FREQ card). The 

format for a typical PARAM card is: 


Col's 1-8 PARAM 

Col's 9-16 one of the 8 computer variable names defined by 
Equation (12) 


Col's 17 


24 real part of variable defined in Col's 9-16 


Col's 25 


32 imaginary part of variable defined in Col's 9-16 (only 
non zero entry is for variable FI0MEG) 


Comments 

Strictly speaking, the real part of variable FI0MEG should be 0.0; however, 
for the NUSC Univac 1108, operating with the level 15.5 version of NASTRAN used 
to implement the procedure, an arbitrary small number if entered (say 1.0x10 ) 

in order to avoid a strange system type error message that is printed when 
exactly 0.0 is entered as the real part. The FI^MEG variable is only used to 
compute and subsequently print out the internal forces at the cuts after all 
the main calculations for displacement are completed. The mentioned error mes- 
sage probably will not appear if other NASTRAN versions and/or other computer 
systems are used. 


DMI Cards for Matrix Partitioning 


A set of DMI direct matrix input cards are needed to provide the informa- 
tion NASTRAN needs to partition the mass, damping and stiffness matrices. 

Three groups of cards are needed; a column partitioning vector for the left cut 
group of displacement node components, CV100; a column partitioning vector for 
the interior group of displacement node components, CV010; and a column parti- 
tioning vector for the right cut group of displacement node components, CV001. 

A set of row partitioning vectors are automatically generated by the Appendix A 
DMAP instructions. The column partitioning are made up of entries that are 
either 1.0 or 0.0. Since all entries within NASTRAN are assumed to be zero 
unless otherwise specified, the user need only enter 1.0 values in the appro- 
priate slot in each of the above mentioned partitioning vectors. The rules are 
simple and are as follows: 

• Formation of left cut partitioning vector CV100 

Enter a 1.0 in each row number corresponding to each active indepen- 
dent component degrees-of-freedom lying along the left cut. The 
length of the CV100 vector is L+I+R and there should be L 1.0 entries 
(the remaining I+R entries are automatically zero by virtue of not 
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being defined). If the left cut nodal numbering pattern is sequential 
and starts with the lowest node number of the whole system (e.g., node 
1. 2, 3, . . .), then the first L entries of CV100 will be all 1.0 
values. However, if the left cut numbering scheme does not contain 
only the lowest node numbers, but instead the whole system is num- 
bered at random, then the L 1.0 entries will correspondingly be dis- 
tributed throughout the CV100 vector, and the "bookkeeping" involved 
with defining the CV100 vector becomes messy. The user having MPC's, 
SPC s or OMIT's applied to nodes along the left cut must be sure to 
account for these during the process of entering the 1.0 values into 
the partitioning vector CV100. 

Formation of the interior partitioning vector CV010 

Enter a 1.0 in each row number corresponding to each active indepen- 
dent component degree-of- freedom lying on the interior of the struc- 
ture. The length of the CV010 vector is L+I+R and there should be I 
1.0 entries (the remaining L+R entries are automatically zero). If 
the interior nodes are numbered sequentially, (following the same 
sequential pattern used in the CV100 vector), then the middle L+l, 

L+2, . . . L+I entries of the CV010 vector will all be 1.0 values. 
Again remember to account for SPC's, MPC's and OMIT's in the number- 
ing scheme. 

Formation of the right cut partitioning vector CV001 


Enter a 1.0 in each row number corresponding to each active indepen- 
dent component degree-of-freedom lying .on the right cut of the peri- 
odic structure. The length of the CV001 vector is L+I+R and there 
should be R 1.0 entries (the remaining L+I entries are automatically 
zero). If the left cut, interior, and right cut nodes are all num- 
bered sequentially (in the respective order mentioned), then the end- 
ing L+I+l, L+I+2, . . . L+I+R entries of the CV001 vector will all be 
1.0 values. Again remember to account for SPC's, MPC's and OMITS's 
in the numbering scheme. 


DMI Cards Format 


The bulk data cards for the definitions of the partitioning vectors via 
the standard DMI cards is as follows: 

• CV100 vector cards 

-Col 1 s, _E.nir.y 

first header card 

1-8 DMI 

9-16 CV001 

17-24 0 
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Col 's Entry 

25-32 2 

33-40 1 

41-48 1 

49-56 blank 

57-64 integer value equal to magnitude of (L+R+I) 

65-72 1 

second card 

1-8 DMI 

9-16 CV001 

17-24 1 

25-32 enter integer row number ,say Nl, of first 1.0 entry 

33-40 1.0 for entry N1 

41-48 1.0 or 0.0 for entry Nl+1 

49-56 1.0 or 0.0 for entry N 1+2 

57-64 1.0 or 0.0 for entry Nl+3 

65-72 1.0 or 0.0 for entry N 1+4 

73-80 Continuation card name, say CONTI, if needed. 

third continuation card (if needed) 

1-8 +0NT1 

9-16 1.0 or 0.0 for entry N1+ 5 

17-24 

25-32 

33-40 

41-48 

49-56 

57-64 

65-72 1.0 or 0.0 for entry Nl+1 2 

73-80 Continuation card name, say C0NT2, if needed. 

fourth continuation card (if needed) 
similar to third continuation card ... etc. 


• CV010 Vector Cards 

These cards are made up analogously to the CV100 vector already described 
above, the only differences being that on the first two cards, replace CV100 
with CV010; also, the number of the first 1.0 entry slot (N1 value in col's 
25-32 of the second CV100 data card) is different. 


• CV001 Vector Cards 

These cards are made up analogously to the CV100 vector already described 
above, the only differences being that on the first two cards, replace CV100 
with CV001 ; also the number of the first 1.0 entry slot (N1 value in col's 
25-32 of the second CV100 data card) is different. 
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Finally, examples of making the partitioning cards is given later in a 
demonstration problem. 


DMI Cards for Merge Operations 


A set of DMI direct matrix input cards are needed to define certain dummy 
matrices which NASTRAN needs to successfully merge certain internal 
matrices within the Appendix-A DMAP sequence. These input cards will consist 
of a group of eight cards that must always be included for a run. The only 
thing that changes from one NASTRAN run to another is the lengths of these 
dummy arrays. These dummy null arrays are only introduced to avoid DMAP error 
message printouts for Uni vac 1108, level 15.5 NASTRAN that occurred when the 
lead matrix entry in the standard DMAP MERGE operation is not defined. 


DMI Cards Format 


Col 's 


1-8 

9-16 

17-24 

25-32 

33-40 

41-48 

49-56 

57-64 

65-72 


Entry 

first header card 
DMI 

LIXLI enter characters, (does not mean multiply LI times LI) 

2 

1 

1 

blank 

enter integer equal to L plus I 
enter integer equal to L plus I 


1-8 

9-16 

17-24 

25-32 

33-40 


second card 

DMI 

LIXLI 

1 

1 

0.0 


Col's 9-16 
Col 's 57-64 
Col 's 65-72 


Col's 9-16 


third card 

same as first card, except 
enter LIXL2 

enter integer equal to L plus I 
enter integer equal to L plus L 

fourth card 

same as second card, except 
enter LIXL2 
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fifth card 

same as first card, except 

Col's 9-16 enter L2XLI 

Col's 57-64 enter integer equal to L plus L 

Col's 65-72 enter integer equal to L plus I 

sixth card 

same as second card, except 
Col's 9-16 enter L2XLI 

seventh card 

same as first card, except 

Col's 9-16 enter L2XL2 

Col's 57-64 enter integer equal to L plus L 

Col's 65-72 enter integer equal to L plus L 

eighth card 

same as second card, except 
Col's 9-16 enter L2XL2 


DMI Cards for Unit Matrix Definition 

Finally, a set of DMI direct matrix input cards as needed to define the 
unit matrix [I u ] employed in Equation (11). 


DMI Cards Format 

Col's Entry 

first header card 

1-8 DMI 

9-16 UMATR 

17-24 0 

25-32 6 

33-40 1 

41-48 2 

49-56 blank 

57-64 enter integer equal to L plus L 

64-72 enter integer equal to L plus L 
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Col's 


Entry 

second card 

DMI 
UMATR 
1 
1 

1.0 

third card 

same as second card, except 

Col's 17-24 enter 2 

Col 1 s 25-32 enter 2 

fourth card 

same as second card, except 

Col 's 17-24 enter 3 

Col 's 25-32 enter 3 


etc. 


last (L+l s ^) card 

same as second card, except 

Col's 17-24 enter integer value equal to L 

Col's 25-32 enter integer value equal to L 


DEMONSTRATION PROBLEM 

The use of the DMAP sequence in conjunction with the new PARAM and DMI 
input cards defined in the previous section will perhaps be better understood 
by including a specific example solution that has the features that (1) is a 
small size problem convenient for matrix operation checking and debugging 
purposes; (2) contains most of the main ingredients typical of a representa- 
tive problem, thus the system has mass, stiffness, and damping; (3) the exact 
solution to the problem is known for checking purposes. The problem illus- 
trated in figure 4a meets all of these conditions, and corresponds to a plane 
pressure wave propagating in an infinite acoustic fluid medium. Upon sampling 
the pressure response along any two parallel vertical cuts (separated by the 
horizontal distance Lp) it can be shown that. the response (pressure and dis- 
placement) is different only by the amount e^P* where for this particular 
problem 


y* = (|)(L p ) Cos* (13) 


1-8 

9-16 

17-24 

25-32 

33-40 
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where w = driving frequency (rad/sec), L p = length (in.) between two parallel 

cuts, c = compressional wave speed (in/sec); ip = angle of incident wave. The 
solution for the pressure and motion response is sought, for the shaded region 
shown in figure 4a, Only the dotted outlined region is modeled with finite 
elements and is correspondingly shown in figure 4b, 

In order to exercise this demonstration problem to the fullest extent, we 
will use different types of boundary conditions on all four sides of the fig- 
ure 4b finite element model. 


Boundary Conditions 
Left and right vertical cuts 

Since this problem falls within the class of problems solvable by 
the phase difference type boundary condition, boundary conditions 
specified by Equations (7) are enforced. Here the left and right 
cut forces and displacements are taken as unknowns. 

Top face 

The boundary condition for the top face is different from the left 
and right face in that here we explicitly apply the free field pres- 
sure (converted into equivalent nodal forces). The formula for the 
freely propagating pressure wave is given by the expression 


P(x,y,t) = p Q e 


i(kr) iwt 


(14) 


P(x,y) 


where k = w/c = wave number (in. - *) 


r = spatial coordinate normal to direction of wave propagation (in,) 
p = steady state pressure amplitude (psi) 


o 

i = 


V 


■1 


p = pressure (psi) (spatial variation) 
and r is related to the x,y coordinates by the relation 
r = x cosip + y sinip 

The y-direction force at upper face node 3 is computed by substituting 
p(x,H) from Equation (14) into the expression 


(15) 


F*) 


‘3'lb 


y 




x=L p /4 


p(x,H) dx 


(16) 
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the force at upper face node 4 is computed similarly by 

x=3L p /4 

( 4 f“) ■ f- p(x,H)dx 

y J 

x-Lp/4 

and finally the upper face node 9 is computed similarly by 


(17) 


x=L 

(■* = 

9 r y 

x=3L / 4 


/■ 


p(x,H)dx 


(18) 


The demonstration problem is evaluated for the following specific input data: 


p Q = 100, psi 

c = 60000o in/ sec 

L = 2o0 in 0 
P 

H = 2c0 in. (19) 


f = w/2rr = 3000c Hz 
ip = 45° 

4 

p = mass density = O 000096 lb.- sec/in 


Substituting the above Equation (19) input constant into Equations (16), 
(17) and (18) results in 


(,F?) = -49.8972 e 1 

y 

( 4 F?) = -99.7944 e 1 


28.6378° 

38.1837° 


( 20 ) 



= -49.9742 e 1 


47.7297° 


An important point must be made regarding loading the final force array 
{F} (Equation (4)) for the periodic structure problem. Observation of the 
Equation (10) loading vector of the new periodic structure problem statement, 

reveals that the left cut applied forces, {F^}, are applied, as normal, to 

a* q 

the corresponding left cut node; also, the interior applied forces, {F^}, are 
applied, as normal, to the corresponding interior nodes; however, the right 
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cut applied forces, (F a ), are not applied to the right cut nodes, but rather, 

are first multiplied by the complex constant e" ^ 1 , and then applied to the 
corresponding left cut node. For the demonstration problem at hand, substi- 
tuting Equation (19) into Equation (13) implies that y* = 25,45584°, thus 

e 1y * ( q F a J = -49,9742 e i22 ‘ 27387 ° (21 

r y 

therefore, in summary, at node number 3 in the y direction apply a net sum 


force = -49,8972 e 1 28,6378 .499743 e i22„27387° 
and at node number 4 in the y direction apply a net 


( 22 ) 


force = -99,7944 ^ 38,1837 (33) 

and at node number 9 in the y direction apply a net 
force = 0,0 (i.e., do not apply any force). 

• Bottom face 

The boundary condition for the bottom face could have been selected 
similar to the top face (i.e„, we convert the pressure into equiva- 
lent nodal forces). However, instead we use a slightly more compli- 
cated boundary condition that permits us to introduce a [B] matrix 
entry into the problem. More specifically, the pressure and normal 
velocity along the bottom cut can be obtained from the expression 

p(x,o) = pcV p (34) 

where is the particlevelocity normal to the wavefront propagation direc- 
tion. Therefore the relation between the V y velocity component and the y 
direction resisting reaction is given by 

(F a ) « ( p( x ,0) *AA) = pcAA Simp V y (25) 

damping 

constantEC . 

d 

Thus, the bottom face boundary condition (simulating an approximate wave 
absorbing boundary condition) is achieved by placing viscous dampers along 
the bottom cut (see figure 4b), wherein the damper constants are defined by 

C d = pcAA Simp (36) 

where AA is an appropriate area factor relating the pressure and concentrated 
force (F..) . When the wave length of the incident wave is long relative to 
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the mesh size, one can set AA = the surface element length for bottom surface 
nodes off the cuts and set AA = 1/2 the surface element length for the nodes 
lying on the cuts. Thus for the demonstration problem at hand, = 8.14587 

for the middle damper and half that amount for the end cut dampers. 


Preparation of Demo NASTRAN Input Data 

• executive control 

The form of the DMAP instructions presented in Appendix A are general and 
are not problem dependent. The only times the user deviates from the presented 
DMAP sequence is (1) when he switches from one level of NASTRAN to another 
wherein the ALTER statement numbers change; or (2) when he wishes to turn off 
or turn on the intermediate matrix printout switch ISW, defined in the third 
DMAP ALTER statement (ISW = +1 prints all intermediate matrix operation steps 
in addition to a printout of the FORVEC vector which lists the left cut, 
interior, and right cut nodal forces; ISW = -1 no intermediate printout). 

• case control 

The standard CASE CONTROL deck is shown in the APPENDIX A, wherein the 
only thing worth noting is the fact that a DLOAD card is used for the purpose 

of superimposing the two vectors {F®} and e~ y 1 {F^} which are both applied to 

the same left cut nodes (in correspondence with the first partition of Equa- 
tion (10)). 


• bulk data 

The CDAMP2 cards are used to define the damping constants applied at the 
bottom surface nodes. The CQDMEM membrane elements (and corresponding MAT2, 
PQDMEM cards) are used to define the fluid media, employing the displacement 
fluid element approach described in ref „ 8 . The collection of DLOAD, DAREA, 
DPHASE, RL0AD1, TABLED1 cards are used to insert the applied forces defined by 
Equations (22) and (23). The 6RDSET card is employed to eliminate the non- 
applicable 3, 4, 5, 6 degrees-of-freedom for the 2-D membrane elements. Stand- 
ard GRID cards define node coordinates and the standard FREQ card defines fre- 
quency of f = w/2tt = 3000 Hz. The special nine PARAM cards defined by Equations 
(12) are evaluated using the data in Equations (19). In the case of the DMI 
cards for the partitioning vectors CV100, CV010, CV001, the first thing to 
establish for the problem is the sizes L, R, I for the individual partitions. 

It is convenient, for bookkeeping purposes, to number the left cut sequen- 
tially; and to also number the interior sequentially (with the lowest interior 
node number appearing next to the highest left cut node number) and finally 
number the right cut sequentially with the lowest right cut node appearing next 
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to the highest internal node number. Thus for the problem at hand^ L = 2 
degrees-of-freedom/node times 3 left cut nodes = 6; I = 2 degrees-of-freedom/ 
node times 3 interior nodes = 6; and R e [_ = 6 due to the periodicity of the 

t !F*7*I h ?». f1rs U ?o!u tries . 0f CV10 ° are U0 > the rest of entries are 
zero, the 7 th through 12th entries of CV010 are 1.0, the rest being zero; and 

finally the 13th through 18th entries of CV001 are 1.0, all other entries beinq 
zero. The 8 merging DMI cards LIXLI ••• L2XL2 are defined according to the 
instructions given earlier in the paper and need no further comment. Finally 
the unit matrix DMI cards for the UMATR matrix are defined according to the 
instructions given earlier in the paper. 


Results of NASTRAN Demo Problem 

Selected output results for the figure 4b demonstration problem are pre- 

irAPPFNnTy A R TRA ?h dl ^ Ct / rin J 0U ^ ir ea L and 1 n a g inary part type output format) 
in APPENDIX B. The first part of the NASTRAN printout illustrates the net 

3^ - ( lr 9 ^ re l ult J 0f Ec < uations ( 22 ') and (23 )); if is always good to 
include as a checking feature of the input. 

Next the stress output is printed and refers to the stress computed at the 
° f the element. Note that due to the sign convention difference regard - 
ing stress and pressure (opposite in sign, see ref. 8 for details), the user 
must reverse the sign of stress to obtain the pressure, i.e.. 


-a 


xx 


- -a 


yy 


/ . P s(l ° u1cl be noted that the NASTRAN output format headings are in error 

ni^cToLt 01 ^* 1119 bu9 in NASTRAN that has remained in practically all levels 
of NASTRAN); the output heading should be read as follows If nr earh element m. 
of stress output): Real Pt. (<,„/, Real Pt. (o yy ). Real Pt. (a xy ), Img.Tt. 

Img. Pt. (<Jyy ) , I mg. Pt. (a xy ). The results shown in APPENDIX B, Table B-i, 

fll2 W -ii he - NASTRAN r es “ 1ts n ext to the exact solution, after converting the real 
and imaginary parts into amplitude and phase data. 


(a xx } ’ 


r fhat h ?n^to\ ?l Ut 1S evalu ^ e f Wlth P(x>y) from Equation (14) at distances 
mTH™Li° C * te 3 l ine i dra Tu P arallel to the wave front) which passes through the 

solut nn T nt * I he N ?? TRAN res ults agree quite well with the exact 

solution in both phase and amplitude. It is noted that a very sliqht asvmmetrv 
exists between the NASTRAN phase angle results for element 3 Y and element 4. Y 
J? robably dua t0 the fact that the boundary conditions on thetopand 

thp h Ln U fnr UrfaCeS d0 + 0t resu1t in exact1y the same applied nodal forces (e.g., 
1 no ? com P uted usln 9 a consistent pressure-to-force applica- 

Ilnt In? 13 Jhat involves the displacement shape functions). A similar com- 
ment applies to the relation between the bottom surface forces and left cut 


^Suppose, for illustrative purposes, node 5 had the x displacement SPC con- 
Str M^ ed $°J er ? and suppose the x displacement of node 4 was forced fthrouoh 
case P - re a i ° n t0 j the - displacement of node 6 - In this particular 9 
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forces. When the program is run with the detailed print switch triggered 
(ISW = +1), the left and right cut forces are computed from the solution dis- 
placement vector and subsequently printed in condensed format (under the FORVEC 
header). Inspection of these results showed that corresponding left Q and right 
cut internal nodal forces differed by the proper amount, y* = 25.455 . 


CONCLUDING REMARKS 

The DMAP sequence presented in APPENDIX A, permits the NASTRAN user to 
solve a class steady state dynamically loaded periodic structures, wherein the 
internal forces and resulting response at the periodic ends of the typical 

i u* 

repeating substructure are related by a known complex phase relation, e 
The DMAP sequence is general and need not be changed from one problem to the 
next. When using the methodology presented here, one should be extremely care- 
ful that the CV100, CV010, CV001 partitioning vectors are prepared properly. 

For large problems, a preprocessor should be written which automatically gen- 
erates these vectors. One should also be careful to allow for SPC, MPC and 
OMIT cards which compact the solution vector and consequently should be taken 
into account in the preparation of the partitioning vectors. The user is 
strongly advised to study the sample demonstration problem presented here 
before undertaking anymore complicated problems. 

At this point, the DMAP sequence has been checked out for a set of small 
size problems, therefore comments regarding computer run time on large problems 
cannot be made at this time. Some of the DMAP operations can be increased in 
efficiency by employing the matrix operation ADD5 DMAP module in place of 
repeatedly employing the matrix ADD DMAP module. The ADD5 routine was not 
used in the 1108, level 15.5 version of NASTRAN due to the fact that the system 
did not consistently successfully add a string of 5 matrices during certain 
checkout phases of the programming. 

Finally, it is advised that before attempting to exercise the DMAP sequence 
presented in APPENDIX A, one should run the demonstration problem as a bench 
mark check in. If a Univac 1108, level 15.5 version is used, it may be necessary 
to increase the "max-files" to 40 in NTRAN$. 
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APPENDIX-A 

(NASTRAN DEMONSTRATION PROBLEM INPUT) 


NASTRAN EXE CUTIVE CONTROL DECK 


ID NUSC PERIODIC STRUCTURE 5AWPLT 

APP PISP 

SOL 8» 0 

T IME 30 

DIAG 2,3»8»14,15 

DIAG 2 2 

ALTER 138 


S DMPA INSTRUCTIONS FOR PERIOnTC FTRUCt UPPC (rongn By A . J- KALINOVSKIY 

$ ISW=+1 PRINT ... =-l NO PRINT OF ALL KEY MATRICES FOR PERIODIC STR. 

PARAM //C.NtNOP/V.N.ISWs+l $ 

PARAM //C»NrNOP/V»N*IMl=-l S 

PARAM //C»N>NQP/V,N.PURSWl=-l g 

$ DEFINE ROW PARTITIONING VECTORS WITH DUMMY ADD OPERATION 

g WHERE CV100 CVOlft CVOOl ARE RFAO IN ON DMT BULK DATA 

ADD CV100./RV100 S 

ADD CV010./RV010 g 

ADD CVOOl »/RV001 g 

COND MAR»ISW g 

MATPRN RV100»RV010>RV001» »// g 

MATPRN CV100»CV0l0*CV001>»// g 

LABEL MAR $ 

g CONDITIONAL PRINT OF KDDrMDP. BDD BEFORE PARTITION APPLIED ( 

COND KAL»ISW g 

MATPRN KDD>MDD»BDDn// » _ 

T" SAVE ORIGINAL KDO MDD BDD MATRICES FOR PROCESSING FORCES AT CUTS 

ADD KDPt/SAVKDD g ' 

ADD MDD./SAVMDD g 

ADD BDP./SAVBOD g 

LABEL KALg 

g PAR IT ION K»M»B MATRICES 

g PARTITION LL BLOCK 

PARTN KDOtCV10O.RVln0/»»»DKLL/C»N,l g 

PARTN MDD»CV100.RV100/»»»DMLL/C*N,1 g 

PARTN BDD»CVlOO»RVlnO/»»>DBLL/C>N.l g 

g PARTITION LI BLOCK 

PARTN KDOtCVOlO»RVlOO/t mOKLI/C»N, 1 g 

PARTN MDO *CV01O*RVlC0/»f» DMLI/C * N. 1 S 

PARTN BD0»CV010»RVl00/>»»DBLI/C»N.l $ 

g PARTITION LR BLOCK 

PARTN KDD»CV001 > RV100/» >> DKLR/C > N»1 g 

PARTN MDD»CVQ01.RV100/m»DMlR/C»N»1 S 

PARTN BDD»CV00l>RVl00/» t »DBLR/C>N>1 g 

g PARTITION IL BLOCK 

PARTN BDD*CV100.RV010/»»tDBlL/C»N»l g 
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PARTN KDDrCVlc0.fTv0ro/7“»»DKiL/C»N*l * 

PART N MDD > C V 100,R VO 10/ > » »DMlL/C » N» 1 $ 

$ PARTITION II BLOCK 

PARTN KDD»CV010.RV010/f » fDKH/C»N»l $ 

PARTN MDD»CV0l0#RV010/<-~fDMlI/C»Nrl $ 

PARTN BDQ»CV010*Rv010/» » fDBlI/C*Nf 1 $ 

$ PARTITION IR BLOCK 

PARTN KDO»CV00l#RV010/f » fDKlR/C»N»l $ 

ParTn moo»cvooi.rvoio/» > >dmir/c»n>i $ 

P ARTN BP O »CV001,Rv01Q/ > , > DBIR/C > N»1 $ 

$ PARTITION RL BLOCK 

PARTN K0D» CV100* RV0Ql/» r »DKRL/C »N» 1 $ 

PARTN M00» CVlOOf RV00i/» t »OMRL/C*N»l $ 

PARTN BOO> CV10C* RVOQl/» • »PBRL/C»N» 1 $ 

S PARTITION Rl RLOCK 

PARTN KCJ»CV010.RV001/o»DKrI/C»N>1 s 

PARTN MODfCVOlO.RvOOl/r * »0MRI/C»N»1 S 

PARTN BDD»CV010,RV001/» , *DBRI/C»N»1 S 

% PARTITION RR BLOCK 

PARTN KDD> CVO01* RV0Ql/» » >DKRR/C»N» .1 $ 

PARTN MDDt CV001# RVOOl/ » t »OMRR/C »N» 1 S 

PARTN BOO* CVOOl* RVOQl/> mOBRR/C > N» 1 $ 

SCONDITIONAL PRINT OF PARTITIONS 

coND kalin »isw $ 

MATPRN DKLL»DMLL»DBLL. .// $ 

matprn dkli»dmli»obli> »// $ 

PiATPRN DKLR»0MLR»D8LR» .// * ~ 

MATPRN DKILrPMIL>DBlL> ,// $ 

MATPRN DKII#OMII»DBlI» »// $ 

MATPRN DKIR>PMIR»DBiR» »// $ 

MATPRN OKRL » DMRL # DBRL » » // $ 

MA TPRN DKRI>PMRI»DBRl» »// $ 

MATPRN DKRR»DMRR*DBRR» *// * 

LABEL KALIN $ 

$ * 

$ FORM PARTITIONS OF ASSEMBLED MATRIX MDOX» BOOX> KDDX 

$ FORM LL BLOCK 
”S 

$ FORM AMLL 

ADD DMll 1 6MLK/D1 AMLl/V"» Y »"P0NE/V » Y » CMST * 

ADD D1aMLL»DMRL/D2AmLL/V»Y»P0NE/V*Y»CMST $ 

ADD D2 AMLL. DMRR/ AMLl/V,Y» PONE/ V»Y. PONE 4 
PURGE KDD*MDD#BDD/PURSW1 % 

PURGE D1AMLL»D2AMLL/PURSW1 S 
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"5 FORM ABLL 

ADO DBLL r PB LR/D1BLL/ V r Y r PONE/V , Y r CMST t 

ADD DlBLL * DMLR/D2BLL/V * Y r PONE/V r Y , NSMSTO $ 

ADD D 2 BLL , 0 K LR/D3BLL/V , Y r PONE/V » Y » SMSBO g 

ADD D3BLL • DBRR/D4BLL/V , Y r PONE/V » Y , PONE S 

PURGE D1ABLL>D2ABLL.D3ABLL/PUR£W1 $ 

ADD D4BLL,0MRL/P53LL/V,Yr PONE/V, Y,SMSTO « 

A0D_D 5BLLi_QKR L/ABl-L/ V j Y r.PQNE/V . Y r NSMSBO S 

PURGE D4BLL r D5BLL/PURSW1 $ 

J FORM AKL1 

ADD DKLLrOKLR/DlKLL/V,Y,PONE/VrY,CMST g 

ADD _ D1KLL»DBLR/D2KUL/V rY r PONE/V, Y , NSMSTO g 

ADD D2KLL,0KRL/D3KLL/V,YrP0NE/VrY,CMST $ 

ADD D3KLL , QKRR/D4KLL/V >Y» PONE/V , Y . PONF $ 

ADD 04KLL r DBRL/AKLL/ V r Y » PQNE/V » Y , SMSTO * 

PU RGE DlKLLrP2KLL.D3KLL,D4KLL/PURSWl g 

$ FORM LI BLOCK 

$ 

* FORM AML I 

ADD DMLIrDMRI/AMLl/VrY, PONE/V, YrCMST * 

S FORM ABLI 

ADD DBLI , D3RI/01BLI/V rYr PONE/V, Y,CMST g 

ADD D1BLI ,DMRI /D2BL I /V,Y» PONE/V, Y, SMSTO * 

ADD D2BLI , OKRI/ABLI/V , Yr PONE/V , Y, NSMSBO g 

PURGE D1BLI»D2BLI/PURSW1 $ 

$ FORM AKLI 

ADD DKLI»DKRI/D1KLI/V»Y, PONE/V, Y, CMST * 

ADD DlKLl r DBRI/AKLI/V # Y t PQNE/V, Y, SMSTO % 

PURGE D1KLI/PURSW1 $ ~~ 

* FORM IL BLOCK 
$ 

$ FORM AMIL 

ADD DMIL*DMIR/AMIL/V*Y, PONE/V, Y*CMST $ 

* EPM A£IJ. 

ADD DBIL»D8IR/D1BIL/V»Y, PONE/V, YrCMST * 

ADD DIB I L,0MlR/D2BlL/V,V» PONE/V ,Y , NSMSTO g 

ADD D2BIL,0KIR/ABIL/V rYr PONE/V , Yr SMSBO S 

PURGE DlBILrD2BIL/PURSWl $ 

$ FORM AKIL 

ADD DKIL tDKIR/DIKIL/V rYr PONE/V r YrCMST S 

XDO t)lKIL,DBlR/AKlL/V , Y r PONE/V »Y, NSMSTO $ 

PURGE D1KIL/PURSW1 $ 

1 Form ii block 

* AMI I SAME AS DMH ABII SAME AS DMlI AKII SAMf AS DKII 

* FORM RL BLOCK 
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ADD UMATK * / ABKL/V7Y * SM5B0 5 

AOD UMaTR*/AKRL/V,Y,CMST $ 

T'FOKM rr block 

ADD UMaTR./AKRR/V,Y»NONE S 

1 THE LR » IR .RI BLOCKS ARE NULL 

$ print blocks before assembly 

COND DE8»ISW $ 

MATPRN AMLL»ABLL»AKLL» >// S 

matprn amli»abli*akli,»// S 

MATPRN AMIL»ABIL»AKIL>»// $ 

MATPRN aBRL»AKRL*AKRR»>// * 

PRTPARM //c»n»os 
LABEL DEB* 

$ ♦ 

$ NEXT ASSEMBLE BLOCKS WITH MERGE 
$ * 

$ MERGE LL BLOCK 

MERGE LIXLI » r » AKLL»CV100»RV100/KDDXLL/C » N* 1 * 

MERGE LIXLI» » *AMLL*CV100rRV100/MDDXLL/C»N»l * 

MERGE LIXLI» » »abll»cvioo»rvioo/bddxll/c»n*i * 

S MERGE LI BLOCK 

MERGE LIXL2» > » AKLI »CV010>RV100/KODXLI/C»N» 1 $ 

MERGE LIXL2* » »AMLI.CV010»RV10n/MDOXLI/C*N.l * 

MERGE LIXL2»»»ABLI»CV010»RV100/8DDXLI/C»N.1 S 

1 FORM PARTIAL SUM 

ADD KDDXLL*KDDXLI/SUMK1 * 

ADD MDDXLL f MDDXLI/SLJMM1 * 

ADD BDPXLL > BDDXLI/SUMB1 $ 

$ PURGE COMPONENTS OF PARTIAL SUM NOT NEEDED ANY MORE 

purge kddxll»kddxli.mddxll>mddxli»bddxll.bddxli/purswi S 

* MERGE il block 

MERGE L2XLI>»»AKIL»CV100»RV010/KDDXIL/C>N»1 * 

MERGE L2XLI»»»AMIL»CV100»RV0l0/MDDXIL/C»N»l * 

MERGE L2XLI m> ABIL»CV100 »RV01O/BDDXIL/C >M» 1 $ 

* CONTINUE PARTIAL SUM 

ADD SUMKl »KDDXIL/SUMK2 $ 

AOD SUMM1»MDDXIL/SUMM2 S 

ADD SUMB1»BDDXIL/SUMB2 < 

COND JOHN» ISW $ 

MATPRN SUMK1 » KDDX IL . SUMK2 , SUMM1 # MDDX IL//* 

MATPRN SUMM2 » SUMB1 p BDDXlL » SUMB2 i //* 

LABEL JOHN $ 

"5 PURGE COMPONENTS OF PARTIAL SUM NOT NEEDED ANY MORE 
PURGE SUMK 1 # KDDX IL » SUMM1 1 MDDX IL » SUMB1 * BDDX I L/PURSW1 s 
MERGE II BLOCK 
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MERGE. L2XL2#**DKII*CV010#RV01n/KDDXII/C»N.l 5 

MERGE L2XL2 > *»DMII>CV010»RV0lQ/MDPXII/C>N»l $ 

MERGE L2XL2»*»DBlI»CV010»RV010/BODXII/C»N.l * 

S CONTINUE PARTIAL SUM 

ADD SUMK2» KDDXI I/SUMK3 $ " “ " 

ADD SUMM2»MDDXII/MDDX $ 

* WHERE MDDX IS THE FINAL MASS MATRIX SUM 

ADD SUMB2 » BDDX 1 1 /SUMB3 S 

COND BULL.ISW * 

MATPRN <^UMK2t KDnXT I « SUMK5» SUMmP » MDDXT T //< 

MATPRN MDDX t SUMB2 » BDDX 1 1 # SUMB3 , //$ 

. L AB E L-BULLS 

5 PURGE COMPONENTS OF PARTIAL SUM NOT NEEDED ANY MORE 

PURGE SUMK2>KDDXlI*SUMM2fMDDXlI,SUMB2.BDPXII/PURSWl « 

* MERGE RL BLOCK 

MERGE 1 TXLI > t » AKRL »rVlQO»RVOQl /KDnXRi /riM.l $ 

MERGE LIXLIm»ABRL»CV100»RV001/BDDXRL/C»N»1 S 

$ CONTINUE PARTIAL SUM 

ADD SUMK3 » KD0XRL/SUMK4 $ 

ADD SUMBStBDDXRL/BDDX S 

$ WHERE BDDX IS THE FINAL MATRIX SUM ' 

COND PILL* ISW $ 

MATPRN SUMK3 » KDDXRL . SUMK4 , SUMB3 » BDDXRL//S 

MATPRN BDDXr »>,//$ 

LABEL PILLS 

S PURGE COMPONENTS OF PARTIAL SUM NOT NEEDED ANY More 

PURGE SUMK3» KDDXRL »SUMB3»BDDXRL/PURSW1 S 
S MERGE RR BLOCK 

MERGE LIXLI*»*AKRR»CV001»RV001/KDDXRR/CfN»l S 
S CONTINUE PARTIAL SUM 

ADD SUMK4»KDDXRR/KDDX $ ' 

$ WHERE KDDX IS THE FINAL SUMMED MATRIX 
COND BOOK » ISW~$ 

MATPRN SUMK4 » KDDXRR , KDDX > » //S 

LABEL BOOK $ 

S PURGE COMPONENTS OF PARTlAI SUM NOT NFFDEn ANY MORF 

PURGE SUMKL#KDDXRR/PURSW1 $ 

S ♦ 

s *** NOW ALL THE KDDX BDDX AND MDDX MATRICES ARE FORMED 

COND JACK » ISW * “ 

MATPRN MDDX *BQDX» KDDX »,// $ 

LABEL JACKS 

ALTER 139 >139 

FRRD C ASEXX » USETD , DLT , FRL , GMD , GOD , KDDX , BDDX » MDDX . » D It/UDVF , PSF , PDF » PPF/ 
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c . N »DT5P/C * N . D IRECT/V » N » LUSEfD/ VV nTmpcFI /V , N » S I NGLE/V #N# OMT . , 
V.N>NONCUP/V,N »frqset/c»y»decomopt=i * 


ALTER 140 
COND ABEfISW $ 


S PRINT THE FORCE VECTORS USED BY THE EQUIVALENT PERIODIC STRUCTURE 
MATPRN PDF.PSF#PPF» .// $ 

* PROCESS THE LEFT AND RIGHT PERIODIC CUT FORCES 


ADD SAVMDC»SAVKDD/DSAV/V»Y»NOMEG2/V»Y.P0NE S 


ADD DSAV * SAV8DD/KNET/V t Y t PONE/VrY • FIOMFG * 

PURGE DSAV/PURSW1 $ 

MPYAD KNET »UDVF, /FORVEC/C »Nf 0/C *n»i/c»n»o/c»n»i * 
MATPRN KNET »UDVF»FORVEC » f // $ 

LABEL ABE $ 

enoalter 


CEND 




DISPLACEMENT s ALL 


UatcHMMW 





sorted sulk data echo 


1 1 

♦ .. 2 

.. . _.3L. 

• t . 4 

• * 

5 6. 

.... J. 

• i 8 — 

.. 9 

.. 10 ... . 


101 

4.0?293 

1 

2 






CDAMP2 

106 

8.14587 

6 

2 






CDAMp? 

107 

4.07293 7 

2 






CQDMEM 

1 

1 

2 

5 

4 

3 




cqdmem 

2 

1 

1 

6 

5 

2 




CODMEM 

3 

1 

4 

5 

8 

9 




CQDMEM 

4 

1 

5 

6 

7 

8 




darea 

21 

3 

2 

-49. 

8972 





darea 

21 

4 

2 

-99. 

7944 





□area 

210 

3 

2 

-49. 

9742 





dload 

10 

i. 

i. 

100 

1. 

101 




DMI 

CV001 

0 

2 

1 - 



18 

1 


DMI 

cvooi 

1 

13 

1. 

1# 

1 . 

1. 

1. 

D06S 

+OGS 

1 * 









DMI 

CV010 

0 

2 

1 

1 


18 

1 


DMI 

CV010 

1 

7 

1. 

1# 

1. 

1. 

1 L 

PIGS 

+ IGS 

1. 









DMI 

CV100 

0 

2 

± 

1 


18 

1 


DMI 

evior 

1 

1 

1. 

1. 

1. 

i. 

1. 

CATS 

+ATS 

1# 









DMI 

L I XU I 

0 

2 

1 

1 


12 

12 


DMI 

LIXLI 

1 

1 

.0 






DMI 

LIXL2 

0 

2 

1 

1 


12 

12 


DMI 

LIXL2 

1 

1 

^0 






DMI 

L2XLI 

0 

2 

1 

1 


12 

12 


DMI 


1 

1 

.0 






DMI 

L2XL2 

0 

2 

1 

1 


12 

12 


DMI 

L2XL2 

1 

1 

_ ^0 






DMI 

UMATR 

0 

6 

1 

2 . 


6 

6 


DMI 

wnrsim 

1 

_i. 

i. - 






DMI 

UMATR 

2 

2 

1. 






DMI 

UMATR 

_ 3 

3 

1 . 






DMI 

UMATR 

4 

4 

1. 






DMI 

UMATR 

5 

5 

_U 






DMI 

UMATR 

6 

6 

1. 






DPHASE 

23 

3 

2 

28.6378 





dphAse 

23 

4 

2 

38,18377 





dphase 

230 

3 

2 

22.27387 





FRETS 

SRDSET 

1 

3000. 





3456 




1 


.0 

.0 

.0 





GRID 

2 


m 

HoHul 

.0 






3 


.0 

2.0 

.6 





GRID 

4 


1.0 

2.0 

.0 





GRID 

5 


1.0 

1.0 

•1 





GRID 

6 


1.0 

.0 

.0 





GRID 

7 


2.0 

Mil 

.0 ... 





GRID 

Q 


2.0 

1.0 

.0 





GRID 

9 


_2.S_ 

2.0 

HI 





MAT2 

10 

345600. 

345600 

. .0 

345600. 



.000096 


PARAM 

CMST 

.9029168.0 







PARAM 

FIOMEG 

1.0-20 

18849. 

55 






PARAM 

N0MEG2 

-3.553+8.0 







PARAM 

NONE 

-1*0 

.0 

















PARAM 

NSMSTO 

-8108.830.0 









+1*0 

__.0 







param 

SMSBO 

+2.280- 








PARAM 

SMSTO 

+8108.830 . 0 







pqomem 

1 

10 

1.0 







RIOAni 


21 



22 





RLOAD1 

101 

210 


230 

22 





TABLED X 

22 








PULS1 

+ULS1 

ENDDATA 

0.0 

1.0 

looooo 

. 1.0 

ENDT 
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APPENDIX B (NASTRAN DEM O NSTRATION PROBLEM OUTPUT) 


r REQuENCY = 3.000000+03 

x_a.ji£_L. 

EX 1 0 4 D 

— M-Ei-g.JLo R 




(REAl/IMASINARY) 




-'UXNT J U. 

TYPfc 

T1 

T2 

T3 

R1 

R2 

R3 


3 

G 

*0 

-9.003838+01 

.0 

*0 

.0 

.0 

. 


. 

.0 _ 

-4+J28562l+ftl_ ... 

.n 

■— t-Q 

.0 

.0 


4 

G 

.0 

-7.844159+01 

.0 

.0 

.0 

.0 




.0 

-6.169147+01 

.r 

.0 

.0 

.0 



-REQUENCY s 3.00ft000+03 - 

COMPLEX ST RE S S E S — I..N Q U AJL BILA TERAL M F MBRaNES (C OD MEM) 

’"THEALMMMpiAffY) — 


element 

ID. 


1 

2 

3 


4 


rea1 °xx ^TORMAt-<it^. rea1 5 _ 

-9.069812+01/* -9.069812+01 

-9.808248+01/* _-9. 808248+01 
-7 . 895863+01/ . -7.895863+01 

-9.080797+01/* -9.080797+01 


stresses - iR~ ELEMENT "coordinate s~ystem“- 

- JLg— q xy TWRMf teaLjBglgjg 1ma 9 Qyy ’tRETm-XXlmag 

•0 /» -4.316528+01 -4.316528+01/* .0 

• 3 “2 1? 073P2+0 1 -2.207309+ 0 1 /, ,0 

•0 /* -6.209007+01 -6.209007+01/, .0 

*® ./ * -4j, 31424 3+01 _ -4.314243+ 01 /, .0 


?jcy 


TABLE B-l 

NASTRAN-EXACT SOLUTION COMPARISON 


p 

element^^ 

NASTRAN 

|P|.Psi 

EXACT 
|P|» Psi 

NASTRAN 
Phase of p 

EXACT 

Phase of p 

1 

100*53 

100.00 

12.683° 

12.728° 

2 

100.44 

100 0 00 

25.451° 

25.456° 

3 

100 * 44 

100.00 

38.180° 

38.183° 

4 

100.53 

100.00 

25.412° 

25.456° 
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b) Infinite Ribbed Plate 



~ H L r 
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F.fg. 3 Cyclic Periodic Structure 







Fluid Media 



Fig. 4 Demonstration Problem 



